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ABSTRACT 



Context. A crucial issue in star formation is to understand the physical mechanism by which mass is accreted onto and ejected by 
a young star. To derive key constraints on the launching point of the jets and on the geometry of the winds, the visible spectro- 
polarimeter VEGA installed on the CHARA optical array can be an efficient means of probing the structure and the kinematics of the 
hot circumstellar gas at sub-AU. 

fH Aims. For the first time, we observed the Herbig Ae star AB Aur in the Ho- emission line, using the VEGA low spectral resolution 

O . (R=1700) on two baselines of the array. 

I Methods. We computed and calibrated the spectral visibilities of AB Aur between 610 nm and 700 nm in spectral bands of 

O 20.4 nm. To simultaneously reproduce the line profile and the inferred visibility around Ha, we used a 1-D radiative transfer code 

Y^ (RAMIDUS/PROFILER) that calculates level populations for hydrogen atoms in a spherical geometry and that produces synthetic 

C/2 spectro-interferometric observables. 

^ Results. We clearly resolved AB Aur in the Ha line and in a part of the continuum, even at the smallest baseline of 34 m. The small 

P-Cygni absorption feature is indicative of an outflow but could not be explained by a spherical stellar wind model. Instead, it favors 

, a magneto-centrifugal X-disk or disk- wind geometry. The fit of the spectral visibilities from 610 to 700 nm could not be accounted 

for by a wind alone, so another component inducing a visibility modulation around Ha needed to be considered. We thus considered 

a brightness asymmetry possibly caused by large-scale nebulosity or by the known spiral structures. 

Conclusions. Thanks to the unique capabilities of VEGA, we managed to simultaneously record for the first time a spectrum at a 
resolution of 1700 and spectral visibilities in the visible range on a target as faint as my = 7.1. It was possible to rule out a spherical 
geometry for the wind of AB Aur and provide realistic solutions to account for the Ha emission compatible with magneto-centrifugal 
acceleration. It was difficult, however, to determine the exact morphology of the wind because of the surrounding asymmetric ne- 
l/~i bulosity. The study illustrates the advantages of optical interferometry and motivates observations of other bright young stars in the 

^^ same way to shed light on the accretion/ejection processes. 

o 

I Key words. Methods: observational Techniques: high angular resolution - Techniques: interferometric - Stars: individual (AB Aur) 

. . Stars: emission-line - Stars: circumstellar matter 
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1. Introduction tain. However, determining if similar processes are at work in 

d these objects is extremely important, as it would help to esta- 

The class of Herbig Ae/Be (HAEBE) stars gathers the pre-main- ^j^^j^ ^j^^^j^^^ ^j^^ accretion/ejection mechanism is universal ir- 

sequence stars of intermediate mass (1.5 < M/Mq < 10), with respective of the mass of the central body The hot gas that can 

spectral types of B to F8, strong emission lines, and the pres- ^^ ^^^^j^^^ -^ accretion and ejection flows close to the source 

ence of infrared to submillimeter excess flux. Like their solar- ^^^ ^^„ ^^ ^^^^ ^^ ^^^^ ^1^^ corresponding physical conditions 

type analogs, the T Tauri objects, HAeBe objects are known to ^j^^^^ j^ -^^ ^^^i^^i^^ -^ ^-^^^^ especially those in the optical range 

be surrounded by protoplanetary disks of gas and dust, responsi- ^j^^^ j^^^^ ^ ^-^ diagnostic potential. These phenomena, how- 

ble for this excess emission. In the case of lower mass stars, the ^^^^^ ^^^^^ -^ ^ ^^^^jj -^^ ^f ^ f^^ astronomical units (AU) 

formation is commonly explained by the gravitational collapse ^^^^^^ ^^^ ^^^^ corresponding, for close star formation regions, 

of a dusty disk cloud and magnetically controlled accretion via ^^ ^ ^^^ milliseconds of arc. Such a smaU scale can only be re- 

a disk. The formation of more massive stars is still very uncer- solved with interferometrv 
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Guided by these considerations, we have attempted for 
the first time to observe a bright visible Herbig Ae/Be star, 
AB Aurigae, with a spectro-interferometer operating at opti- 
cal wavelengths, the VEGA spectrograph (Mourard et al. 2009) 
installed at the CHARA Array (.Brumme laar et al. 2005 ). This 
unique combination of spectral and spatially resolved informa- 
tion allowed us to study the physical conditions for the gas 
producing the Ha line on sub-AU scales. Interpreting the re- 
sults, however, required appropriate modeling of the emitted ra- 
diation. Using a radiative transfer code, we obtained observa- 
tional insights on the generation mechanism of winds around 
AB Aurigae, giving experimental support to the widely ac- 
cepted, but not yet tested, theory of magneto-centrifugal acce- 
leration. In Sect. 2, we recall the present knowledge about the 
target. In Sect. 3, we describe the observations and the data pro- 
cessing. In Sect. 4, we present the radiative transfer code and 
discuss the wind models adopted to interpret the observed Ha 
visibilities and spectrum. 

2. The target: AB Aurigae 

AB Aurigae (AB Aur; AOVe; d=144 pc) is the brightest 
Herbig Ae star in the northern hemisphere, and is often con- 
sidered as the prototype of the HAeBe class. Its photometric and 
spectral properties are well studied over a wide wavelength range 
"l5hm & Cat ala 1993 Catalaetal. 19991 [Grady et al. I999a 



Table 1. Journal of AB Aur observations on October 8 and 9, 
2008. 



IVan den Anc ker et al. 1999T l and there is evidence of both 
accretion and outflows in AB Aur JCatala & Kunasz 19871 
Garcia-Lopez et al. 2006 IWassel 20091 1. The Ha line exhibits a 



P-Cygni profile, which is a common probe of stellar wind, and 
is known to be variable, especially in the blue wing of the line. 
Among the alternatives evoked to explain the observed variabi- 
lity were a wind coming either from the equatorial regions of 
the star or from a circumstellar rotating disk dCatalaet al. 1999l l. 
The wind in AB Aur seems to be complex, probably controlled 
by a magnetic field, even if no magnetic detection of AB Aur has 
been reported so far jAlecian et al. 2007T l. 
AB Aur is surrounded by a complex combination of gas- 
rich and dust-dominated structures probed on large scales by 
imaging. Using the HST, Grady et al. (1999b) detected ex- 
tended nebulosity and circumstellar material at all position an- 
gles out to 1320 AU in the R-band, whereas Fukagawa et 
al. (2004) detected a double spiral structure at 200-450 AU 
in the H-band. Spatially resolved observations at high angu- 
lar resolution at near-infrared, thermal infrared, and millimet- 
ric wavelengths have shown the presence of a circumstellar disk 
jMillan-Gabet et al. 20011 Eisner et al. 20041) . rotating at non- 
Keplerian velocities ( Pietu et al. 2005I I. All these observations 
agree on a low inclination for the system (i.e., i<40°), in dis- 
agreement with the value of 70° + 20° deduced from spectro- 
scopic studies JCatala et al. 19991 ). A strong near-infrared con- 
tinuum emission was found inside the silicate dust sublimation 
radius (0.24 AU) and was interpreted as due to a hot gaseous 
disk jTannirkulam et al. 200 8 ). A localized asymmetry on spa- 
tial scales of 1-4 AU has also been detected by Millan-Gabet et 
al. (2006) who speculated that it might result from either loca- 
lized viscous heating due to a gravitational instability in the disk 
or from a close stellar companion, in agreement with spectro- 
astrometric observations in the optical lines jBaines et al. 2006l l. 

3. Observations and data processing 

Data were collected at the CHARA Array with the VEGA spec- 
trometer recording spectrally dispersed fringes at visible wave- 
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B(m) 
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2008-10-08 
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33.9 
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34.0 
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-lh22 
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HD 29646 


33.9 
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-4.94 
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2008-10-09 


08:21 


HD 29646 


86.1 


120.9 


-2h32 


2008-10-09 


08:57 


AB Aur 


89.9 


119.4 


-2h32 


2008-10-09 


09:25 


HD 29646 


95.8 


112.0 


-lh43 


2008-10-09 


09:49 


AB Aur 


99.0 


110.3 


-lh25 



lengths. Two pairs of telescopes were combined: W1W2 (whose 
projected baseline B~100 m) and S1S2 (B~34 m). Observations 
were performed around the Ha line, between 610 and 700 nm, at 
the lowest spectral resolution of VEGA (R=1700). Observations 
of AB Aur were sandwiched between those of a nearby calibra- 
tion star (HD 29646) that was chosen to be bright enough (my 
- 5.7) and to have a similar spectral type (A2). The observation 
log is given in Table [T] 

Each set of data was composed of observations following a 
sequence calibrator-star-calibrator, with 10 files of 3000 short 
exposures per observation. Each data set was processed using 
the Ci estimator and the VEGA data reduction pipeline de- 
tailed in Mourard et al. (2009), with a large spectral band of 
20.4 nm (shifted by steps of 10.2 nm over the observing range) to 
reach a sufficient signal-to-noise ratio. The resulting resolution 
of the spectral visibilities is R=30, leading to an interferometric 
field of view of 240 mas. The poor signal-to-noise ratio meant 
that the visibilities could not be processed for the data recorded 
withWlW2. Only the squared visibilities coiTesponding to S1S2 
were computed in 9 spectral channels and calibrated using an an- 
gular diameter of 0.20±0.01 mas for the calibrator HD 2964^ 
(Fig.[T}right). 

The VEGA spectrum of AB Aur exhibits an Ha emission 
line with a clear P-Cygni profile (Fig. fll-left). The line to con- 
tinuum ratio is at most about 2.5, the line FWHM is 6 A, and 
its emission peak is slightly red-shifted. In addition, the line 
presents a very broad absorption from ~640 nm until ~665 nm, 
probably due to Stark broadening in the absorption of the radia- 
tion by nearly hydrostatic cold layers located just above the star's 
photosphere ( Felenbok et al. 19^31 ). Due to the width of the pro- 
cessing spectral window, the squared visibilities at 646 nm, 
656 nm, and 666 nm include all or part of the Ha emission, 
and the values of the squared visibility are between 0.72 and 
0.75. Although they are relative to continuum emission alone, 
the squared visibilities at 636 nm and 676 nm are below unity, 
while those farther away from Ha are consistent with 1. AB Am- 
is thus clearly resolved by the S1S2 baseline both in the Ha Une 
and in a part of the visible continuum. 

4. l\/lodeling and discussion 

The spectroscopic P-Cygni feature in AB Aur shows that the 
line is formed at least in part in a wind. Therefore, we focused 
on a wind model to constrain the properties of the emitting re- 
gion. To simultaneously reproduce both the Une profile and the 
spectral visibilities, we have adapted two numerical codes in 
spherical symmetry for the production and transport of hydro- 
gen (H) lines ( [Rajabi 2010) . The first one, RAMIDUS, calcu- 
lates the level populations of H in a moving fluid for assigned 
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Fig. 1. Left. Ha line recorded at R = 1700 with VEGA. Right. Calibrated squared visibility at R = 30 obtained with the S1S2 
CHARA baseline. 
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Fig. 2. Best model predictions compared to the VEGA spectrum (left) and spectral visibilities (right), both in full circles with error 
bars. Dashed lines correspond to the X-wind model and full lines to the disk-wind one. For the squared visibilities, triangles denote 
the model prediction including both the wind and the contribution of the asymmetric reflection nebula. 



radial dependencies of the basic wind parameters (density, ve- 
locity, temperature). The code uses the Sobolev approximation 
and the escape probability method JCastor 1 9701 1 . In cascade, the 
PROFILER code calculates synthetic images, line profiles and 
visibilities of H lines for various geometrical configurations. By 
selecting limited portions of the spherical winds, we approxi- 
mated 2-D axisymmetric models. We did not include accretion 
in the models. 

A full spherical wind providing matter along the line of sight 
can be readily excluded, because it would produce too much 
blue-shifted absorption with respect to the observed P-Cygni 
feature that is shallow but not very deep (normalized flux>0.88). 
Even if the issue of the inclination angle of the AB Aur system is 
still alive (see Sect. 2), we chose a wind geometry in agreement 
with a small inclination (i< 40°), and thus excluded a double 
cone geometry, reminiscent of a bipolar jet, which would lead, 
for small inclinations, to a line profile with two distinct emis- 
sion features instead of the observed single peak. Inspired by the 
magneto-centrifugal (MC) scenario for the acceleration of stel- 
lar jets, we assume that, on the observed scales (much less than 
2-3 AU from the star/disk), the wind occupies a region delimited 
by a flattened torus defined by the pale gray region in Fig. [3] 



This configuration can simulate the base of either a X-wind 
(ICamenzi nd 1990 Shu et al. 1995) or a disk-wind (hereafter D- 
wind; IFerreira 199 7 , Konigl & Pudritz 2000| l accelerated along 
the magnetic field lines anchored in the inner few AU of the ac- 
cretion disk. In this region the wind has a wide angle aperture, as 
it has not yet reached the Alfven surface, after which magnetic 
collimation occurs and the wind collimates into a jet. In Fig. [3] 
R,=0.02 AU is the corotation radius, Rj=0.24 AU is the inner 
radius of the optically thick dusty disk, R/=2.3 AU is the exten- 
sion of the wind, and a is the angle between the disk mid-plane 
and the surface of the wind, constrained to be < 60° for MC ac- 



celeration to occur ( Spruit 1996) . In our X-wind model, all the 
material in the flow comes from the region located very close 
to the corotation radius and has a constant mass outflow rate 
M„, through the torus section, while for the D-wind M„, progres- 
sively increases from R, to about 1 AU, to simulate that the wind 
originates in an extended region of the disk. The wind velocity 
is assumed to be V(r) = Vq + Vf{l - (R*/rf), where Vo = 20 
km.s"' is the initial velocity, r the radial distance from the star, 
R^, the stellar radius, and Vy the terminal velocity. The density 
profile is dictated by V{r) and M„(r) through the mass continuity 
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Table 2. Parameters of the models, with Nn.ej the average hydro- 
gen density at the wind footpoints. 



system, already reported by several authors ([Grady et al. 1999b| 



Model 



Vf [km/s] M„ [Mp/yr] NH.ej [cm-^] 



X-wind 0° 
D-wind 20° 
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Fig. 3. Geometry adopted for the empirical X-wind and disk- 
wind models described in Tablel2] R,, R^;, and Rj are the corota- 
tion radius, dust sublimation radius, and extension of the wind, 
respectively. 



equation, while the gas temperature profile T(r), which reaches 
a maximum of 12 000 K, is chosen in accordance with the the- 
ory of wind heating illustrated by Shang et al. (1998) and Garcia 
et al. (2001) for the X-wind and D-wind, respectively. Modeling 
of the Stark absorption in the extended line wings has been in- 
cluded using a template spectrum from Kurucz (1979). 

Regarding the modeling of the line profile, both an X-wind 
and a D-wind could provide a good fit (Fig.|2[left) for the model 
parameters collected in Table l2] To reproduce at the same time 
the line width and the absence of a strong blue-shifted absorp- 
tion, the best predictions for the X-wind have nearly a null incli- 
nation (/) and a large opening angle (a) of 60°. For the D-wind 
model, since the disk hides part of the high velocity material 
in the red-shifted wing, the best predictions are obtained with 
a small inclination (20°) and a smaller opening angle of 35° to 
avoid too much blue-shifted absorption. Both models lead to a 
small inclination in agreement with the previous interferometric 
determinations. As pointed out by Pogodin (1992), a change in 
the opening angle of the wind (a) affects the depth of the blue- 
shifted absorption. Concerning the mass loss rate M„ , our best 
prediction is of the same order as previous determinations from 
line profile analysis JCatala & Kunasz 19871) for the X-wind but 
is significantly higher for the D-wind. In fact, our various sim- 
ulations clearly show that the estimate of M„ is very dependent 
on the velocity law and the temperature profile considered in the 
wind model, and on the scale on which the wind is accelerated. 
Since the geometry of our D-wind model strongly deviates from 
a spherical wind, the value of Mw obtained for this model cannot 
be directly compared to the previous determinations of Mn. 

Regarding the modeling of V^, the X-wind provides only a 
negligible decrease in Ha (V^ = 0.98) with respect to the con- 
tinuum (Fig. |2]-right). In this case, the emission is concentrated 
at R,, while it is much more extended in the D-wind model that 
produces a larger decrease in V^ (down to 0.89). In both cases, 
however, the observed visibility variation could not be modeled 
with the wind alone so we considered an additional component 
producing an extended continuum emission. This component 
could come from the extended reflection nebulosity around the 



Fukagawa et al. 2004 IMillan-Gabet et al. 20061 ). In practice, we 
used this secondary source to reproduce the modulation of the 
visibility in the continuum. Our wind models combined with this 
additional emission provided a good fit of the visibilities. For in- 
stance. Fig. l2 right corresponds to a secondary source separated 
from the star by 38.2 mas (i.e. 5.5 AU) contributing to 7% of the 
total flux. Because of the large number of free parameters and the 
fact that we only have measurements along one baseline direc- 
tion, we found many different configurations for the secondary 
emission that could fit, for both wind models, the observed visi- 
bilities (e.g. a much larger separation of 591.2 mas). 

In conclusion, thanks to the unique capabilities of 
VEGA/CHARA, we managed for the first time to record si- 
multaneously a spectrum at a resolution of 1700 and a visibi- 
lity spectrum in the visible range on a target as faint as my = 
7.1. Despite the limited signal-to-noise ratio of our data, they 
allowed us to rule out a spherical wind and to make realistic 
predictions on the nature of the Ha emission region at sub-AU, 
and our data are compatible with a magneto-centrifugal mech- 
anism for the production of the wind. It was difficult, however, 
to determine the exact morphology of the wind and disentan- 
gle the X-wind and the D-wind because of the extended nebu- 
losity. Moreover, all models neglect the contribution of mass 
accretion to the H-line emission, so self-consistent models of 
winds including this process are needed. Complementary obser- 
vations with larger CHARA baselines are required to place more 
stringent constraints on the complex environment and the photo- 
sphere contribution of AB Aur that can be resolved in the visi- 
ble with VEGA. We also expect to take advantage of the new 
instrumental facilities of CHARA to study other bright HAeBe 
and provide promising insights on the accretion/ejection mecha- 
nisms in intermediate-mass young stars. 
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